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In the literature on integrable systems we find Hamiltonian operators 
without explanation. There is a notable silence on Lagrangians from which 
these Hamiltonian and symplectic structures can be derived. We show that 
starting with Lagrangians, which turn out to be degenerate, the Hamiltonian 
operators for integrable systems can be constructed using Dirac's theory of 
constraints. We illustrate this by giving a systematic discussion of the first 
Hamiltonian structure of KdV. First by Dirac's theory and then applying 
the covariant Witten-Zuckerman theory of symplectic structure we arrive at 
its flux. Then we turn to a new Lagrangian for KdV recently obtained by 
Pavlov and derive the corresponding new symplectic structure for KdV. We 
show that KdV admits infinitely many Lagrangian formulations and therefore 
infinitely many symplectic structures. 

1 Introduction 

While there is no precise definition of complete integrability there are many 
properties that we expect from a completely integrable system which in fact 
enable us to recognize it as such. Bi-Hamiltonian structure for which we have 
the celebrated theorem of Magri is a case in point. Non-linear evolution 
equations that can be cast into Hamiltonian form in two inequivalent but 
compatible ways admit a full set of conserved quantities required by complete 
integrability. These conserved Hamiltonians are in involution with respect 
to Poisson brackets defined by both Hamiltonian structures. 

We recall from earliest school that Hamiltonian structure is derivable from 
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a Lagrangian, yet in complete integrability there is a notable, even deafening 
silence on Lagrangians. Why? 

It turns out that in the variational formulation of integrable equations the 
Lagrangians are degenerate and therefore require the use of Dirac's theory 
of constraints p| in order to achieve Hamiltonian form. On the other hand, 
degenerate Lagrangians for integrable systems always admit a complete set 
of second class constraints which enables us to eliminate all the momenta 
and write the Hamiltonian equations of motion solely in terms of the original 
variables. In effect this means that we can often guess the Dirac bracket with- 
out recourse to Dirac's theory. Thus in the literature on integrable systems 
Dirac's theory of constraints hardly even gets a nod while Dirac brackets 
are generally known as Hamiltonian operators. Even in the excellent book of 
Dorfman that she was able to complete before her untimely death, Dirac's 
name appears in the title but in its pages no reference can be found to a sin- 
gle degenerate Lagrangian for which the Hamiltonian structure is obtained 
through Dirac's theory of constraints. So, for both reasons, namely the fact 
that Dirac's theory is not in the arsenal of most mathematicians as well as 
the convenience afforded by second class constraints of eliminating the mo- 
menta completely, the present literature contains very little on degenerate 
Lagrangians for integrable systems. It was, however, my original motivation 
for going into integrable systems |^ . 

Inevitably something as fundamental as the Lagrangian would make a 
spectacular come-back and indeed it has. The covariant approach to sym- 
plectic structure which has relatively recently been formulated by Witten 
and Zuckerman employs the Lagrangian as its head-piece. 

The usual Hamiltonian approach is not covariant because it singles out an 
independent variable as "time" with respect to which the evolutionary system 
will be defined. Dirac was well aware of this short-coming because he was 
ultimately interested in the Hamiltonian formulation of Einstein's general 
relativity and in his book [Q] complains several times about the fact that 
the Hamiltonian formalism is not covariant. The Witten-Zuckerman covari- 
ant theory of symplectic structure is very much in the spirit of Dirac's work. 
It starts with a Lagrangian and in the construction of the symplectic 2-form 
the most crucial role is played by the boundary terms in the first variation of 
the action. We shall now illustrate the prominence of the Lagrangian in the 
covariant theory of symplectic structure by applying the Witten-Zuckerman 
construction of the symplectic 2-form to the grandmother of all integrable 
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systems, namely KdV. We shall show that not only do we recover some famil- 
iar results but we shall also able to present some new symplectic structures 
starting from new Lagrangians for the KdV equation. Shockingly enough, 
we shall find that there are still unknown fundamental results in the theory 
of KdV. 



2 KdV as bi-Hamiltonian system 

The bi-Hamiltonian structure of KdV was the first one to be discovered and 
has consequently served as the model for the multi-Hamiltonian structure of 
all integrable systems. It is well known that the KdV equation 

Ut + QuUx + Uj:xx = (1) 

can be cast into the form of Hamilton's equations 

Ut = Mu) = {u,H}o = JS^H (2) 

where X is the vector field defining the flow which is manifest from eq.(|l|), 
J is the Hamiltonian operator defining the Poisson bracket and 6u denotes 
6/Su, the variational derivative with respect to u. The Hamiltonian operator 
which is a skew-symmetric matrix of differential operators satisfying the Ja- 
cobi identities is simply obtained from the Dirac bracket as the subscript D 
in eq.(0) indicates. However, this was not the historical route to the Hamilto- 
nian formulation of KdV which started with an important paper by Gardner 




and independently by Zakharov and Fadeev ITT] where they showed that 
the infinite set of conserved Hamiltonians Hi,i = 1, . . . , oo commute with 
respect to Poisson brackets defined by 

= 0, (3) 

which, as we shall soon show explicitly, is of course a particular example of 
the Dirac bracket. Thus 

Jo = ^ (4) 

ax 

is the first Hamiltonian operator for KdV. The next important development 
was an unpublished but widely known result of Lenard, namely the recursion 
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operator for KdV. Magri realized that it leads to the second Hamiltonian 
operator for KdV 

r. d d ^ 
Ji = — + 2u — + —2u 5 

dx-^ dx dx 

and was able to formulate his theorem on bi-Hamiltonian structure. Lenard's 
recursion operator is simply 

n = J^{Jo)-' (6) 

and in the bi-Hamiltonian formulation of KdV the recursion relation for 
conserved Hamiltonians is given by 

Ut„ = Jo SuHn+l = Jl SuHn (7) 

with Hamiltonian densities given by 

-Ho = \u\ (8) 

Til = u'-\u^ (9) 

5 1 
^2 = -u^-huu^+-uj (10) 

which determines the higher flows in the KdV hierarchy. 

In effect, there was no derivation of these Hamiltonian operators. When 
we are presented with a result like this all we can do is to check the properties 
of skew-symmetry and Jacobi identities required of a Poisson bracket. Then 
we must keep quiet, but it is an uneasy quiet. 



3 Dirac bracket is the Hamiltonian operator 

The first correct identification of the Hamiltonian operator as the Dirac 
bracket appeared in a paper by Macfarlane which, as far as I know, is 
unpublished and has not been given the attention that it deserves. The 
first systematic application of Dirac's theory of constraints to the degenerate 
Lagrangian for KdV was given in . 
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We start with the variational formulation of KdV. For this purpose we 
need to introduce 

u = (p^, (11) 
the Clebsch velocity potential. The potential KdV is 

(j)t + 3(j)^ + (j),^^ = (12) 

and it can be directly verified that the equations of motion following from 
the variational principle 

SI = 0, I = j Cdtdx (13) 

with Lagrangian density 

^^o = \m. + ^1-\4>.1 (14) 

yield KdV. This Lagrangian is degenerate because its Hessian 



(15) 



vanishes identically. Alternatively, the canonical momentum given by 



TT 



2 



(16) 



cannot be inverted for the velocity 0t and we have a degenerate Lagrangian 
system. Therefore the Hamiltonian formulation of the Lagrangian (|T^) re- 
quires the use of Dirac's theory of constraints. 

Following Dirac we introduce the primary constraint that results from 
the definition of the momentum (|16D 

<^ = Ti-\ct>x (17) 

and a priori require that it vanish either "weakly," or as a "strong" equation 
depending on whether or not the constraint is first, or second class in Dirac's 
terminology. This is determined by evaluating the Poisson bracket of the 
constraints. If the result should vanish modulo the constraints then it is 
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a first, otherwise second class constraint. A second class constraint is an 
equation that can be solved to eliminate a canonical variable. 

In order to calculate the Poisson bracket of the constraints we use the 
canonical Poisson bracket relations 

{T^{x)A{y)}p = K^-y) (18) 

with all others vanishing. The result 

{<l>(x), ^y)}p = Uy{x ~y)- U^{y - x) (19) 

shows that the constraint (|1^) is second class as it does not vanish by virtue of 
the constraint itself. A word of caution, that I owe to my friend Galvao [|12| , is 
in order here: One should not simplify the Poisson brackets of the constraints 
(|19D using the rules for manipulating distributions at intermediate stages of 
calculation. In our case we shall need to invert this distribution and therefore 
we must keep the cumbersome appearance of eg. (|T^) . Simplifications should 
be made only at final results. 

The total Hamiltonian of Dirac is given by 



En 



J {^(Pt- C + X^)dx (20) 
/ + I <I>J + A (tt - i^.) ] dx (21) 



where A is a Lagrange multiplier. The condition that the constraint is main- 
tained in time 

mx),HT} = (22) 

gives rise to no further constraints which would have been secondary con- 
straints. Instead, using eqs.(|l^), we find that the Lagrange multiplier is 
completely determined from eq.(p2D 



A = -3 (j)^ - (f)xxx (23) 

which is expected, since the constraint and therefore the total Hamiltonian is 
linear in the momenta the correct equation of motion (^2|) will result only if 
the Lagrange multiplier is simply the component of the vector field defining 
the flow for potential KdV. The total Hamiltonian density of Dirac 

Ht = I <P! - TT {3 (P^ + (24) 
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follows from the substitution of the Lagrange multiplier ( p3D in eq.(^). Now 
the check that all the Hamiltonian equations of motion are satisfied with 
the Hamiltonian (53) is straight-forward. We can summarize all of them in 



Hamilton's equations 

At = {A,HT} (25) 

where A is any smooth functional of the canonical variables 0, vr and their 
derivatives. 

There is, however, one further and very important simplification that 
we can carry out because in Dirac's theory second class constraints hold as 
strong equations. This fact enables us to eliminate the momentum in the 
total Hamiltonian (p^) using the solution of eq.(|16]). Thus we find 



= - 0i + 2 <PJ (26) 

for Dirac's total Hamiltonian density. Apart from an overall minus sign, this 
is just the Hamiltonian function (^ for the first Hamiltonian structure of 
KdV. 

The Dirac bracket is the projection of the Poisson bracket from phase 
space onto the hyper-surface defined by the constraint. Given any two dif- 
ferentiable functionals of the canonical variables A and i3, the Dirac bracket 
is defined by 

{A{x),B{y)}o= {A{x),B{y)} 

r (27) 

- / {A{x),^{z)}J{z,w){^{w),B{y)}dzdw 



where J is the inverse of the matrix of Poisson brackets of the constraints. 
The definition of the inverse is simply 

j {$(x), ^z)}J{z, y) dz = 5{x - y) (28) 

which results in a differential equation to be solved for J. Starting with the 
Poisson bracket relation (|TUp we find that eq.(^) can be solved readily to 
yield 

J{x,y) = J{x)6{x-y) (29) 




-1 



6ix -y) = e{x - y) (30) 
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where 6 is the Heaviside unit step function and (d/dx) ^ is the principal 



value integral [13|. With the definition of the Hamiltonian operator J{x) 
given above this principal value integral is the first Hamiltonian operator of 
KdV in terms of the Clebsch potential. 

Now we need to make contact with the first Hamiltonian operator @ for 
KdV which is expressed in terms of the velocity field u rather than its poten- 
tial 0. In view of the definition of the potential ([Tl|) and the transformation 
properties of Hamiltonian operators, we find 

J d f d \ ^ d d 

^^^^ dx \dx J dx dx 

which is the same as (^). This is the derivation of the first Hamiltonian 
operator for KdV from first principles, i.e. from a Lagrangian approach. 
I have gone into it in some detail not only because it is a nice illustrative 
example but also because carrying out the same construction for the second 
Hamiltonian operator (|^) is a highly non-trivial unsolved problem. 



4 Symplectic form of KdV 

The Hamiltonian operator maps differentials of functions into vector fields 
and in the opposite direction we have the symplectic 2-form u which is the 
principal geometrical object in the theory of symplectic structure. The sym- 
plectic 2-form is closed 

Suj = 0, (32) 

which is equivalent to the Jacobi identities satisfied by the Dirac bracket. By 
Poincare's lemma, in a local neighborhood, u can be written as 

uj = 6a, (33) 

where a is a 1-form. The statement of the symplectic structure of the equa- 
tions of motion consists of 

ij^u = 5H (34) 

which is obtained by the contraction of the symplectic 2-form u with the 
vector field X defining the fiow. 
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The symplectic 2-form is given by 



a; = - 1 5<P{x) A K{x, y) 5<P{y) dy (35) 

where K{x, y) is the inverse of J{x, y). This route of arriving at the symplec- 
tic 2-form through the Hamiltonian operator is commonplace in the literature 
of integrable systems. But in the Lagrangian approach it is absurd because 
we have the answer already in eq.(^): Hamiltonian operators are Dirac 
brackets for systems subject to second class primary constraints which are 
obtained by inverting the Poisson bracket of the constraints! Hence 

ir(x,y) = {$(x), $(!/)}, (36) 

the symplectic 2-form density can be obtained directly from the Poisson 
bracket of second class constraints. In the case of the first Hamiltonian 
structure of KdV from eq. ( pi?D we find 

uJQ = 5(i)^ S^, (37) 

for the first symplectic 2-form of KdV. 

It remains to check Hamilton's equations in the symplectic form ( ^If ) for 
KdV using this symplectic 2-form and Hi. For this purpose we need to 
contract the 2-form (|37|) with the vector field 




defining the flow for potential KdV. Recalling 

i5/5<i>{y)Hx = S:r{x - y), (39) 
it follows immediately that eq.(0) is satisfied with Hi given by (||). 

5 Witten-Zuckerman 2-form 

Time plays a privileged role in Hamiltonian mechanics. While this presents 
no problem for systems with finitely many degrees of freedom, in field theory 
it has the severe disadvantage of non-covariance. In order to remedy this 
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situation Witten |^ and Zuckerman have introduced the conserved cur- 
rent 2-form which provides an elegant covariant formulation of symplectic 
structure. The Witten-Zuckerman 2-form u;^, where in our case /z ranges 
over two values t and x, satisfies 

6uJ^' = 0, (40) 
= 0, (41) 

the properties of closure and conservation. Hamilton's equations can now be 
written in the covariant form 



= SH^" (42) 

where if* is the familiar Hamiltonian function and is its flux. 

The Witten-Zuckerman current 2-form u'^ is derived from the variational 
principle underlying the equations of motion. In the case of the Lagrangian 
( p^ we need to consider its first variation assuming the equation of motion 
(|12D and its Jacobi equation 



6(j)t + 6 5(f)^ + 5(pxxx = 0. (43) 
Then the first variation of the Lagrangian reduces to a conservation law 

6C = a\, + a% (44) 



where is a 1-form. For the Lagrangian (|1^) we find 



dC^A,+ dC^._fdC\ (45) 



■'XX 



■'XX ^Yx- 



and the symplectic current 2-form is given by eq . (^31) . We find that 

uj' = i 50, A 50 (46) 
cu'' = 3(f)x5(f)x A 5(f) + 5(f)x A 5(f)xx + -5(l)xxx /\ 5(f) (47) 
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where we have used the Jacobi equation p3|). Note that the freedom of 
adding an arbitrary divergence term to the Lagrangian disappears in the 
Witten-Zuckerman symplectic 2-form. We can check that eq.(^) remains 
invariant under the replacement of —^(pj^ in the Lagrangian ([1^) by its di- 
vergence equivalent ^(j)x(j)xxx- 

It is obvious that both components of uj^ satisfy the closure property (^) 
but the check of conservation law (|4l|) requires use of the Jacobi equation 
(^31) again. Finally, we need to check Hamilton's equations in the covariant 
symplectic form of eqs.(|4^). We have already checked the t component in 
section ^. Using ( pQ] ) and its generalization to higher derivatives, we find 
that the contraction of (|47|) with (^) yields 6Hf where 

9 1 

Hi = - 0i + 3 (p^Cpxxx -6(px (Pxl - + - (pJx (48) 

is the familiar flux of Hi. 

We conclude by remarking again that the time component of the Witten- 
Zuckerman 2-form P6| ) is precisely the symplectic 2-form ( pTf ) obtained from 
Dirac's theory of constraints. 



6 Pavlov's new KdV Lagrangian 

We have derived the first Hamiltonian operator for KdV as the Dirac bracket 
for the degenerate Lagrangian (0) and showed that it can be obtained from 
the covariant Witten-Zuckerman theory as well. But now we know that there 
is also a second Hamiltonian operator for KdV. So the question naturally 
arises as to 'What is the Lagrangian for which Magri's Hamiltonian operator 
is the Dirac bracket?" 

I have asked this question, thinking that it is a rhetorical question, in 
several talks that I gave on integrable systems. Last time I did so, at NEEDS 



2000 in Gokova, my friend Pavlov |jT^ got up to say that there are indeed 



many Lagrangians for KdV and showed me that 

^-1 = {<Px + I 0...) <Pt + [l ^x - 5 0. 0i + I 0.') (49) 

is another Lagrangian for KdV. Indeed one can verify it directly. But the 
point is that this Lagrangian is not just an inspired guess. There is a direct 
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derivation of ( ^91 ) which opens up the flood gates to an infinite family of 
Lagrangians for KdV as the subscript —1 indicates.^ 

The derivation of Pavlov's new Lagrangian for KdV is based the Lenard 
recursion relation (^. To understand this, let us first go back to the classical 
Lagrangian (0) and see how, in retrospect, we would derive it. For n = 1 
the recursion relation (|^) can be written in the form 

ut = {4>x)t = Jo SuHi = 6^ Hi{u = (f)^) (50) 



in view of eqs.(^) and (|TT|). Indeed, in the Lagrangian (|T4|) terms that do 
not depend on the velocity precisely make up Hi. In order to obtain the full 
Lagrangian ( |T^ ) all that remains to be done is to rewrite the left hand side 
of (|50|) , namely (ptx in variational form. In complete analogy, for Pavlov's 
Lagrangian terms which do not involve the velocity consist of H2 given by 
(p!OD, the next conserved Hamiltonian function in the KdV hierarchy. Then 
we go back to the Lenard recursion relation and write it for n = 2 



Jq^uH2 = Ji6uHi = TlJoSuHi 



(51) 



and therefore for Pavlov's Lagrangian the factor in front of (pt is simply 
obtained by writing the action of the recursion operator on (ptx in variational 
form. From this construction it is manifest that starting with Hn all we need 
to do in order to obtain a new Lagrangian £_„ for KdV is to write TZ^ (ptx in 
variational form. Hence we have 

Theorem: For every Hamiltonian function in the KdV hierarchy, there 
exists a degenerate Lagrangian that yields KdV as its Euler- Lagrange equa- 
tion. 

The number of Lagrangians for KdV is therefore infinite in number. How- 
ever, here we must note that while KdV will be an extremum for the first 
variation of all these Lagrangians, the Euler equations require something 
weaker. Namely, the n^^ power of the action of the recursion operator on 
KdV should vanish. 



have since received email from Z. Popowicz who has written a computer program 
to generate Lagrangian after Lagrangian for KdV. Because he is a computer expert he 
couldn't resist the temptation to invent his own notation. Unfortunately, I found his 
notation absolutely impossible to decipher. Fortunately, there is no need to do so. In the 
next paragraph we shall present the derivation of all these Lagrangians. 
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Finally, we note that the Lenard recursion operator is non-local and its 
repeated application on (ptx will require the introduction of non-local terms 
in higher Lagrangians. This is not a problem because the original KdV 
Lagrangian is itself non-local! The Lagrangian (0) is not expressible in 
terms of the original variable u but by its potential ([TT|) , or = {d/dx)~^ u. 
Higher Lagrangians for KdV can be written in local form by introducing 
potentials for the Clebsch potential according to (j) = ip^ etc. 



7 New symplectic structure of KdV 

Every new Lagrangian for KdV will give rise to a new symplectic structure. 
This can be obtained either through Dirac's theory of constraints, or directly 
through the Witten-Zuckerman theory. We shall now briefly discuss both 



of these procedures for Pavlov's Lagrangian (|49|) . The result will be a new 
symplectic structure for KdV. It is not the one we were questing after: The 
new symplectic 2-form is not the inverse of Magri's Hamiltonian operator. 



7.1 Dirac constraint analysis 

Pavlov's Lagrangian (^) is degenerate because it is linear in the velocity (pf 
Its Hessian (^) vanishes identically. So we need to apply Dirac's theory of 
constraints in order to cast it into Hamiltonian form. The constraint obtained 
from the definition of momentum that follows from the Lagrangian (^9]) is 
given by 

= ^ - 0x - 2 '^^^^ (52) 

and it is second class because 

1 1 

$(y)}p = -S^^^{y - x) - -Syyy{x - ?/) + 2 0^ S^{y -x)-2(py 8y{x - y) 

(53) 

does not vanish modulo (0). The symplectic 2-form, obtained through (^) 
and (^) is given by 

uj_i = ^ 6^^^^ A 50 + 2 0^. 50^. A 50 (54) 



which is manifestly a closed 2-form. Now Hamilton's equations (|3^) are 
satisfied with H2 on the right hand side. 
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In order to obtain the Hamiltonian operator that corresponds to this 
new symplectic structure for KdV we should invert the Poisson bracket of 
the constraints ( ^5] ) through the definition (P5|). Here we encounter another 
surprise because the equation satisfied by the two-point function J_i{x,y) 
that will lead to the new Hamiltonian operator through the definition in 



eg. (|29|) is given by 



d 

— J-i{x,y) +4 0^— J_i(a;,y) + 2 (p^^ J_i{x,y) = 6{x - y) (55) 
dx-^ dx 

which is precisely the equation to be solved for the inverse of Magri's second 
Hamiltonian operator (^. The inverse of Magri's operator is not known. 
One could write it as an infinite series which, however, does not appear to 
be tractable. So we are unable to write the minus first Hamiltonian operator 
explicitly but we have arrived at a curious result. Namely, the inverse Magri's 
Hamiltonian operator is also a Hamiltonian operator for KdV. 



7.2 Witten-Zuckerman analysis 

The covariant symplectic 2-form for Pavlov's Lagrangian is derived along the 
same lines as in section |^. Here we shall only record the final result. The 
time component of the Witten-Zuckerman symplectic 2-form is the same as 
the expression in eq.(0) obtained from Dirac's theory. Now, however, we 
also have its flux 

uj'L^ = {12 + 5 (f)xxx) 5(f)x A S(f) - Q(f)xx5(f)xx ^ 5(f) + 5 (j)x6(f)xxx A 6(j) {5Q) 

1 3 1 

"2 ^'P^x A 50 - 4 0^ 5(f)xx A5(j)x- - Hax a 50^. + -5(l)xxx A Scpxx- 

The check that Hamiltonian equations in the covariant form of eqs.(^) are 
satisfied with this symplectic 2-form and the second Hamiltonian function 
H2 (0) is lengthy but straight-forward. For the space component of the 
Witten-Zuckerman 2-form it is even lengthier to check that these equations 
are satisfied with (|56D and 

= 12 0i^ + lO0,30,,.,-45 0,Vi + 8 0,0j, + 5 0,^^,,^ 

+0XXX05X - 2 04x5 (57) 

the flux of H2 on the right hand side of (0). 
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8 Conclusion 



The Lagrangian approach to integrable non-hnear evolution equations is one 
where we can derive everything from first principles. We had earlier used 
it for various equations ranging from dispersive water waves |jl5| to Monge- 
Ampere equations [l^, in differential geometry. Here we have tried 
to illustrate this approach using the first Hamiltonian structure of KdV. 
Starting with the classical Lagrangian for KdV we applied Dirac's theory of 
constraints to arrive at its first Hamiltonian operator. We showed that the 
covariant Witten-Zuckerman theory of symplectic structure not only yields 
the symplectic 2-form obtained from Dirac's theory but also gives us the flux 
component of the symplectic 2-form. Then we applied these approaches to 
Pavlov's new Lagrangian and arrived at a new symplectic structure for KdV. 
We have arrived at the remarkable result that the inverse Magri's Hamilto- 
nian operator is also a Hamiltonian operator for KdV. Pavlov's procedure 
for constructing new Lagrangians is applicable to all conserved Hamiltonians 
of which KdV has infinitely many. Therefore there will be infinitely many 
symplectic structures for KdV. 
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